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ABSTRACT 

We exploit the strong lensing effect to explore the properties of intrinsically faint and compact galax- 
ies at intermediate redshift (z s ~ 0.4 — 0.8) at the highest possible resolution at optical wavelengths. 
Our sample consists of 46 strongly- lensed emission line galaxies discovered by the Sloan Lens ACS 
(SLACS) Survey. The galaxies have been imaged at high resolution with HST in three bands (Vhst, 
2"8i4 and #i6o)> allowing us to infer their size, luminosity, and stellar mass using stellar population 
synthesis models. Lens modeling is performed using a new fast and robust code, klens, which we test 
extensively on real and synthetic non- lensed galaxies, and also on simulated galaxies multiply-imaged 
by SLACS-like galaxy-scale lenses. Our tests show that our measurements of galaxy size, flux, and 
Sersic index are robust and accurate, even for objects intrinsically smaller than the HST point spread 
function. The median magnification is 8.8, with a long tail that extends to magnifications above 40. 
Modeling the SLACS sources reveals a population of galaxies with colors and Sersic indices (median 
n ~ 1) consistent with the galaxies detected with HST in the GEMS and HUDF surveys, but that are 
(typically) ~ 2 magnitudes fainter and ~ 5 times smaller in apparent size than GEMS and ~ 4 mag- 
nitudes brighter than but similar in size to HUDF. The size-stellar mass and size-luminosity relations 
for the SLACS sources are offset to smaller sizes with respect to both comparison samples. The clos- 
est analog are ultracompact emission line galaxies identified by HST grism surveys. The lowest mass 
galaxies in our sample are comparable to the brightest Milky Way satellites in stellar mass (1O T M ) 
and have well-determined half light radii of 0'/05 (« 0.3 kpc). 



1. INTRODUCTION 

Understanding how galaxies were formed and how they 
evolve into those we see today is an important cosmolog- 
ical question. In hierarchical galaxy formation, gas con- 
denses and cools within a dark matter halo. To form disk 
galaxies, tidal torques impart angular momentum to the 
dark matter halo and associated baryons; angular mo- 
ment um is conserved as the d isk galaxy forms within the 
halo (jFahiLEfstathiou 1980). The most massive late- 
type galaxies are predicted to have formed from mergers 
of smaller progenitors. However, such models are com- 
plicated by the details of star formation, feedback pro- 
cesses, cluster interactions, and effects of the bulge (e.g. 
IMo et al.|[l998l ). 

Models for galaxy formation and evolution predict 
certain relations between the basic physical properties 
of galaxies (i.e. luminosity, size and mass); quantify- 
ing these relations can help test the standard paradigm 
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of galaxy formation and place limits on future mod- 
els. The size-magnitude (or luminosi ty) and size-mass 
relations are well- studied locally (e.g. IS hen et all 120031 : 
iDriver et al. 2005). The relations for early and late- type 
galaxies (typically defined as having Sersic index n > 2.5 
o r < 2.5) are f ound to diverge. 

iShen et all ([2003) looked at the mass-size and size- 
magnitude relations for galaxies in the Sloan Digital Sky 
Survey (SDSS). A characteristic mass of M* = 1O 1O 6 M 
delineates two regimes for the SDSS size-mass relation: 
above this mass, the relation is steeper and tighter 
(r e ff(kpc) oc M^' 39 , <JinR — 0.34 dex) than for less mas- 
sive galaxies (with r e ff(kpc) oc M^' 14 , cri n R — 0.47 dex). 
The behavior of the size-magnitude relation is similar. 

Surveys of intermediate (0.1 < z < 1) and high 
[z > 1) redshift galaxies have attempted to extend stud- 
ies of the size-mass and size- magnitude relations (e.g., 
Ferguson et al.ll2004l:lBarden et al.ll2005l: iMdntosh et al.l 
20051 : iTruiillo et all 120061 : iMelbourne et al.l 12007ft . 
GEMS, an HST survey, allowed iBarden et aT] (|2005f ) 
to study the magnitude-size and mass-size relations 
of late- type galaxies out to z ~ 1. They found a 
dimming of ~ 1 mag from this redshift to z = in the 
rest-frame V-band, but noted that the mass-size relation 
stays constan t. Us ing the Hubble Deep Field South, 
ITruiillo et all ([2006h extended the SDSS and GEMS 
work out to z ~ 2.5, for the most luminous and massive 
galaxies. The authors find that for low Sersic indices, 
galaxies at a given luminosity were ~ 3.0 ± 0.5 times 
smaller at z ~ 2.5, while galaxies at a given mass were 
~ 2.0 ± 0.5 times smaller. 

Me lbourne et al.l (|2007l ) inferred the size-magnitude 
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relation for blue gala xies in the Great Ob servatories 
Origins Deep Survey (jGiavalisco et al.ll200i GOODS). 
These authors found ~ 1.6 magnitudes of dimming in 
the B-band since z ~ 1 for large and intermediate-sized 
galaxies (r e ff > 3 kpc), in agreement with GEMS. Small 
galaxies, on the other hand, were found to have dimmed 
significantly more, by some 2.55 ± 0.38 magnitudes in B; 
this significant evolution is hypothesized to be the result 
of the fading of the starburst galaxies rather than strong 
evolution of the entire small galaxy sample. 

These studies are limited by the resolution and com- 
pleteness of the HST surveys: for GEMS, these lim- 
its correspond to galaxies with M* > 1O 1O M and 
My < —20, as deter mined by the h i ghest redshift bin 

- 1). Meanwhile, ITrujillo et all ([20061 ) could only 
look at galaxies with Ly > 3.4 x 10 10 /i _2 L© and stellar 
mass M* > 3 x 10 10 h- 2 M Q . 

One method of reaching to lower luminosities is to 
use exceptionally lon g exposure time s: the Hubble Ultra 
Deep Field fHUDF: iBeckwith et al.ll2006f ) is the prime 
example of this. Much of the work on the HUDF has 
been centered around high redshift galaxies (z > 1), but 
several authors have looked at the low redshift sample. 
This sample is primarily comprised of faint objects, since 
the field was directe d away from nearby bright galax- 
ies. iCoe et all (|2006[ ) observe large numbers of faint blue 
galaxies with magnitudes as low as Mb = — 14 at z = 0.7 
- believed to be young starburst galaxies - which peak 
at z ~ 0.67. Two additional spectral energy distribu- 
tion (SED) templates with a steep rise towards bluer 
wavelengths were added to accommodate these objects. 

1 Cameron fc Driverl (j2QQ7h look at the size-luminosity re- 
lation for HUDF galaxies with 0.2<z<1.15. The best- 
fit evolution scenario is one with a brightening of 0.9 mag 
and a 5% decrease in size from z ~ 0.1 to z = 0.675. This 
is consistent with the 1 mag arcsec -2 of dimming since 
z ~ 1 found in other surveys. 

Gravitational lensing is another method of extending 
surveys to potentially smaller, fainter, and less massive 
galaxies and does not rely on extremely deep imaging. In 
strong lensing, a massive foreground galaxy deflects the 
light from a background object, resulting in multiple im- 
ages of the source being seen. The source, in addition to 
being distorted, is typically magnified by a factor of ~ 10. 
This phenomenon allows the study of objects smaller 
than otherwise possible: the tiny source galaxy of gravi- 



tation al lens J0737+3216 was s tudie d by [M arshall et al 



( 20071) . while iStark et al.l (|2008D andlSwinbank fc others 



( 20091 ). for example, have used lenses to carry out de- 
tailed studies of high redshift galaxies. 

In this paper, we study a sample of gravitationally- 
lensed galaxi es, selected from the Sloan Lens ACS Sur- 
vey fSLACS: lBolton et al.l 120081 : 1 Auger et al.ll2009l . here- 
after Paper V and Paper IX). These objects were pre- 
viously modeled in the F814W filter; here, we perform 
multi-filter modeling, which allows us to reconstruct the 
source galaxy spectral energy distribution (SED), and 
hence infer its stellar mass. Our aim is to investigate the 
size-mass and size- magnitude relations for these galaxies, 
and thus explore the potential of gravitationally lensed 
galaxies to further the study of galaxy formation and 
evolution. 

This work is organized as follows. We introduce our 



lens sample in Section [2j then give our multi-filter source 
and lens models for the SLACS sources in Section 03 We 
discuss the properties of the lensed sources and make a 
comparison to the GEMS and HUDF samples in Sec- 
tion [H In Section 14. 5| we compare our galaxies to those 
found in recent emission line surveys. We discuss and 
summarize our results in Section [5] A discussion of our 
new lens- modeling code, designed to perform fast and 
robust lens modeling on large numbers of images, is re- 
served for the Appendix; we also present tests on simu- 
lated gravitational lenses and on both real and simulated 
no n- lensed galaxies here. 

Throughout, we assume a flat ACDM cosmology with 
h = £T /(100 km s" 1 ) = 0.7 and ^ m = 0.3. Ap- 
parent magnitudes are given in the AB system unless 
otherwise stated, while Johnson V and B-band abso- 
lute magnitudes are given in the Vega system. All sizes 
are effective radii, a nd we follow IPeng et al.l ([2002) and 
iBarden et al.l (|2005f ) and define them on the major axis: 
for an elliptically-symmetric surface brightness distribu- 
tion, the elliptical isophote containing half the total flux 
has major axis r e ff and subtends area Trqr^ where q is 
the ellipse axis ratio. 

2. THE SLACS LENS SAMPLE 

Here we introduce our sample, a set of galaxies at red- 
shift z s ~ 0.4 — 0.8 being multiply-imaged by massive 
galaxies lying at z^ ~ 0.2. 

2.1. Sample selection 

We use a subset of lenses from the SLACS survey 
([Bolton et al.ll2008l ) as our sample. We select those lenses 
which were classified as "definitely a lens" and were im- 
aged in the Vhst, hi4, and H±qq bands. The instruments 
and filters these refer to are as follows: for the Vhst- 
band, ACS/WFC F555W or WFPC2 F606W, for the 
7 8 i4-band, ACS/WFC or WFPC2 F814W, and for the 
#i6o-band, NICMOS/NIC2 F160W. For the majority of 
the observations a full orbit's exposure (t exp ~ 2000 sec) 
is available, except for 15 lenses that have only HST 
Snapshot images in ACS/WFC F814W (t exp = 420 sec). 
For details on the observations and data analysis, as well 
as full object coordinates, see Paper IX. All images were 
drizzled onto a pixel scale of 0'/05. This multi-filter sam- 
ple comprises 46 lensed galaxies. 

2.2. Subtraction of lens galaxy light 

When performing lens modeling, we attempt to fit the 
often-faint lensed images, which can be hidden by or con- 
fused with light from the lens (foreground) galaxy; thus, 
it is necessary to remove the light of the lens galaxy prior 
to modeling. As the SLACS survey preferentially selects 
bright lens galaxies this is particularly important for our 
work. W e use the radial B-sp line technique, first intro- 
duced bv iBolton et al.l (|2006l in the appendix), for this 
purpose. We refer the reader to Paper V for a full discus- 
sion of this method, but provide a brief summary here. 
First, zero- weight pixels, neighboring objects and poten- 
tial source galaxy features are masked. Second, the data 
is fit using only monopole, dipole and quadrupole terms 
in the angular structure. The lensed features are masked 
based on the residual image from this initial fit; the im- 
age is then re-fitted using higher order multipole terms 
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as necessary. iMarshall et al.l (|2QQ7[ ) estimated that the 
systematic errors in the source size and brightness due 
to the subtraction of lens galaxy light by this method 
are approximately 2% (0.01 kpc for a 0.6 kpc source at 
z s =0.6) and 0.10 magnitudes respectively. This is an 
important but usually not dominant source of system- 
atic error. Similarly, subtraction of the lens galaxy may 
introduce a systematic uncertainty on the Sersic index, 
similar to those introduced by improper sky subtrac tion; 
based on the detailed work by Marsha ll et al.l ([20071 ) , we 
estimate this systematic uncertainty to be of order 0.1- 
0.2. This is larger than random errors, but does not effect 
the results of this paper. 

3. THE SLACS SAMPLE: MASS MODELS AND 
MEASUREMENTS OF SOURCE OBSERVABLES 

In the Appendix, we demonstrate our ability to mea- 
sure galaxy sizes and magnitudes through galaxy-scale 
gravitational lenses with KLENS; here we use klens to 
measure the size, magnitude, and Sersic index of the 
SLACS sources themselves. We refer the interested 
reader to Paper V and Paper IX for images of the SLACS 
lenses prior to lens galaxy subtraction and the recon- 
structed source planes in the F814W filter. 

3.1. Models of lens and source galaxies 

We have modeled the SLACS lenses in three filters; 
the data and model lens planes for the primary modeling 
filter are shown in Figure [U In order to obtain robust 
and easy-to-interpret constraints on the source galaxies, 
we make three fairly standard assumptions, for which we 
give both a priori justification and a posteriori valida- 
tion: 

Assumption 1: As in the testing program described 
in the Appendix, we use only one source galaxy, 
even though in a few cases the source is evidently 
more complicated (see iBolton et al.l |2008[ ). Without 
this assumption, we would be unable to interpret the 
size and stellar mass of the galaxy. The alternative 
would be to use a more complex sour ce model, per- 
haps defined on a pixelated g r id (e.g. I Warren fc Dvd 
2003; Treu & Koopmans 2004; B rewer fc Lewis! I20Q6F 
Suvu et all 120061 : iVegetti fc Koopm ans 2009f h and then 
derive the size and magnitude of the source from this 
complex model. However, imposing a simple, one- 
component Sersic model is standard in the study of 
faint galaxies, even though they may have irregular 
morphologies or mu ltiple knots of star formation (e.g. 
Ferguson etafll200l . If we are to compare our sample 
with others in the literature, we need to derive the same 
parameters as were studied by previous authors. A po- 
tential concern is that oversimplification in the source 
surface brightness distribution may lead to small bi ases 
in the mass model. However, as we discuss in Section [321 
comparison of our mass model parameters with those 
inferred using complex multicomponent source distribu- 
tions shows that, at least for the SLACS sources, this is a 
negligible source of error. In conclusion, considering that 
using a simple analytic source model greatly speeds up 
computation time (<~ 15 minutes in our approach com- 
pared with several hours for a typical pixel-based source 
reconstruction, with a standard desktop CPU), we adopt 
this procedure. 



Assumption 2: We require that the mass model be 
identical across the different filters, in order not to bias 
our source reconstruction. Although ideally this would 
be avoided naturally because there is only one deflector 
potential, lens galaxy subtraction and differences in con- 
trast and signal-to-noise ratio (S/N) can result in small 
differences, which we quantify as follows. Using the igu 
band as reference, we found the difference in <jsie to 
be, on average, 0.6 km s _1 with an r.m.s. scatter of 
13.5 km s _1 in the Hiqo band, and —2.0 km s _1 with an 
r.m.s. scatter of 11.5 km s _1 in the Vrst band. Like- 
wise, the difference in mass distribution axis ratio q was, 
on average, —0.007 with an r.m.s. scatter of 0.15 in the 
Hiqq band and —0.03 with an r.m.s. scatter of 0.12 in 
the Vhst band. Two or three outliers (defined as having 
|A<tsie| > 50 km s _1 or |Ag| > 0.6) were excluded from 
each of these calculations. Fixing the mass model to that 
of the highest S/N reconstruction ensures that we do not 
introduce any unnecessary extra scatter in the properties 
of the source. 

Assumption 3: We fix the source morphology parame- 
ters (position angle, inclination, size and Sersic index) to 
the best-fit values from the filter with highest S/N. (For 
two-thirds of galaxies this filter is F814W; for the remain- 
der, only snapshot images were available in F814W and 
the filter with the highest S/N is F606W. Those objects 
with only snapshot images in F814W are indicated in 
Table [TJ) This is necessary because the S/N sometimes 
differs significantly between filters, which can cause low- 
surface brightness features to be missed and result in 
overly small sizes and Sersic indices, which in turn af- 
fects the inferred magnitude. Thi s procedure, analogous 
to the SDSS model magnitudes ( Ab azajian et al.ll2009l 
and references therein), is effectively equivalent to mea- 
suring colors within fixed aperture, and is widely adopted 
in order to obtain robust colors when the S/N varies sig- 
nificantly between filters. 

Inferred, unlensed angular sizes and apparent mag- 
nitudes for the source galaxies are given in Table [TJ 
along with their lens and source redshifts, and a flag 
indicating whether the Vhst band magnitude refers to 
WFPC2/F606W or ACS/F555W. We also give the to- 
tal magnification ji of each system. This quantity varies 
widely, between 2.4 and 44.3; the median magnification 
is 8.8^5*1, where the error bars correspond to the 16 and 
84 th percentiles. The lens galaxy properties - mass axis 
ratio, inclination, and the velocity dispersion - are not 
presented here, as our concern is with the source galaxies. 
For a thorough treatment of the SLACS lens galaxies, we 
refer the reader to Paper V and Paper IX. 

3.2. Testing mass models: comparison to previous work 

The F814W im ages of the lenses in our sample have 
been modeled by IBolton et al.l ([20 08) and I Auger et al.l 
(|2009h using multiple objects in the source plane. This 
gives us the opportunity to test klens' ability to model 
real lenses and estimate systematic (modeling) errors. 
Although we cannot compare source models (given the 
different definitions), we are able to look at the effects of 
using only one source object on inferred mass model pa- 
rameters. The lens mass model is described by the mass 
axis ratio (g), the inclination, the velocity dispersion of 
the best-fitting singular isothermal ellipsoid (<tsie) and 
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Notes: AB magnitudes are from the Sersic model fits in the source plane, and so are unlensed. The effective radius r e ff is the major axis 
of the elliptical isophote containing half the total flux. Uncertainties on r e ff and the photometry are 13% and 0.3 magnitudes, respectively. 
For full coordinates, see Paper IX. A star in the instrument column indicates that only a snapshot was available in lsi4(F814W). 



the mass centroid. We compare mass axis ratios and ve- 
locity dispersions and find that they agree to within levels 
expected from systematic errors: the average offset be- 
tween our fits and those from Paper Vis— 0.015 ± 0.002 
for q and 2.0 ± 0.2 km s _1 for <jsie- The comparison for 
velocity dispersion is shown in Figure [2j 

We can use the scatter of these relations to estimate 
total errors, including systematics: we calculate 0.07 for 
the error on q and 1.78% for the error on <jsie- These 
are slightly larger than the errors adopted in Paper V 
(0.05 on q and 1.0% on <tsie), as expect ed because our 
model s are less flexible. As discussed by Mar shall et al.l 
(|2007l ) and in the summary at the end of the Appendix, 



if the potential is not perfectly describe d by a SIE (e.g. 
iKoopmans et al.l20Q9t lAuger et al.ll2010f ) then additional 
errors will apply to the source galaxy properties. These 
are taken into account in the analysis of the source pop- 
ulation. 

We note that one of the systems, J0737+3216, 
has a published source s ize, magnitude and stellar 
mass ([Marshall et al.l |2007| ) . Even with a completely 
different lens modeling code, we find lens parameters, 
source size and source Sersic index for this system that 
are consistent with those in this previous work, providing 
further confirmation of the robustness of our approach. 
However, we infer a source that is 1.5 magnitudes fainter 
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Fig. 1. — Lens-subtracted HST images with klens model-predicted arcs. In each case the image in the band with the highest signal/noise 
ratio is shown - this is the image that was used when fitting the lens model. Images are 6 x 6" . 
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uniformly in all bands. T his was traced t o an error in 
the modeling code used bv iMarshall et all (|2007[ h which 
summed the flux in twice-subsampled pixels rather than 
averaging it: this factor of 4 in flux translates to a magni- 
tude difference of 1.5. After correcting for this "bug" the 
magnitudes are also in excellent agreement. The result s 
here supercede those published in Marshall et al.l (|2QQ7[ ). 

These comparisons show that our models are reliable 
and give answers consistent with those obtained by other 
methods once systematic errors are taken into account. 

3.3. Rest-frame luminosity and stellar mass of the 
source galaxies 

We now use our multi-band photometry of the SLACS 
sources to infer rest frame luminosity in the B and V 
(Johnson Vega) bands as well as stellar masses (M*). 
This analysis is based on our models of the unlensed 
source galaxies, where modeling has been carried out 
in three separate bands. Fo r this purpose we u se the 
Bayesian code developed by I Auger et al.l (|2QQ9f ) to fit 
stellar population synthesis (SPS) models to our multi- 
band source galaxy photometry. The code computes, for 
each galaxy, the likelihood of SPS models as a function 
of stellar mass, age, metallicity, star formation history 
and dust content. In combination with a prior on each 
of these parameters, the likelihood then gives a poste- 
rior PDF for each model parameter. The same models 
and posterior can also be used to generate self-consistent 
rest-frame lumi nosities. For this applicat ion we adopt 
SPS models by iBruzual fc Charlotl (|2Q03l ) and uniform 
priors on the logarithm of stellar mass, age, metallicity, 
and time scale of the exponential star formation history, 
as is appropriate for the situation where we do not know 
the order of magnit ude of these quantit ies. 

As discussed by I Auger et al.l ([20091 and references 
therein), although parameters such as age and metallic- 
ity are often degenerate, stellar masses and luminosities 
in the range of wavelength covered by the data can be 
derived quite accurately for a given IMF. Typical errors 
on the transformations to rest frame luminosities are of 
order 0.05-0.1 mags, while typical errors on stellar mass 
are of order 0.1-0.2 dex (see Table [2] for details). For sim- 
plicity we neglect the impact of emission lines on broad 
band photometry, which is estimated to be of order a 
few percent for the typical Ha fluxes (the strongest line 
in the wavelength range of interest) of a few 10 -16 erg 
s _1 cm -2 . Typically, no strong emission line is present 
in the i^i6o-band filter, which is providing the bulk of the 
information for stellar mass estimate. The main residual 
source of uncertainty is the normalizatio n of the initia l 
mass function. In this paper we adopt the lKroupal ([20011 ) 
normalization of the initial mass function (IMF), to fa- 
cilita te comparisons with the GEMS work ( Bar den et al.l 
l2QQ5h . The uncertainty in the lens ma ss density pro- 
file sl ope is an additional source of error ([Marshall et al.l 
120071 ). 

As a sanity check we compared our estimated stellar 
masses wi th those inferred by applying the "standard" 
recipe bv iBell et all (j2003l ) to our inferred colors. For 
the same IMF the stellar masses agree very well, with an 
average offset of 0.02 ±0.02 dex (r.m.s. scatter 0.15 dex). 
Our synthetic rest frame photometry agrees well with 
that inferred by standard "K-correction" procedures (for 
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Notes: Absolute magnitudes are given in the Vega system. Stellar 
masses were inferred assuming a Kroupa IMF. The typical uncer- 
tainty on the absolute magnitudes is 0.17 mag, (0.1 statistical, 0.14 
systematic) with an additional 0.26 mag systematic error (which 
does not affect the color) due to the uncertain lens mass density 
profile slope. This is also the source of a systematic error of 0.1 dex 
on the stellar mass, which is included in the given error bars. 
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Fig. 2. — Comparison of best-fit klens velocity dispersions to the 
Paper V or Paper IX value. SLACS values are on the x-axis and 
the difference between klens and SLACS is shown on the vertical 
axis. The average difference is 2.0 ± 0.2 km s — 1 , and using the 



r.m.s. scatter we estimate the error on (Jsie to be 1.78%. 
error bars shown are total errors (2.5%) on the difference. 



The 



example, by comparing with the V band magnitudes in- 
ferred via a Sbc template we find an average offset of 
0.06 ± 0.06 magnitudes). We conclude that our stellar 
mass and synthetic photometry are robust and unbiased 
within the errors. 

4. PROPERTIES OF THE SOURCE POPULATION 

In this Section we study the properties of the source 
population. To put our galaxy population in context, we 
use two comparison samples of non-lensed galaxies ob- 
served with HST at comparable redshift, by the GEMS 
and HUDF collaborations. After a brief description of 
the two samples (Section 14. ip . we investigate the distri- 
bution of observed properties, i.e. Sersic index, magni- 
tude, and effective radius ( Sect ion |4~3|) . In Section l4~4l we 
investigate the bivariate distribution of rest frame quan- 
tities, i.e. V-band magnitude, B-V color, effective radius, 
and stellar mass. Finally, in Section l4~5l we compare the 
properties of the SLACS sources to those of emission line 
galaxies selected by blind spectroscopic surveys. 

4.1. GEMS-selected comparison sample 

For the fi rst comparison sample, we follow 
iBarden et al.l (|2005[ ) and select, from the publicly- 
available GEMS catalog|3 galaxies with successful 
GALFIT structural fits (GEMS_FLAG = 4) matched within 
0'/5 of a COMBO- 17 object with successful photometric 
redshift estimate (C0MB0_FLAG = 3). We did not select 
galaxies based on their size or Sersic index (other 
than to reject objects with large size or Sersic index 
uncertainties, or that reached the GALFIT boundary 
conditions, as did lBarden et al.ll2005f ). but we did reject 
objects not detected in both F606W and F850LP filters. 

9 The GEMS catalog used in this work can be obtained from 
http://mpia.de/GEMS/gems_20090526.fits, or from M. Barden on 
request. 



The resulting GEMS sample comprised 6999 galaxies 
with measured absolute magnitudes and sizes. We 
adopt the F850LP-measured sizes, since these differ 
from the rest-frame V-band sizes by typically only 3% 
(jBarden et al.l l2QQ5h . We then select galaxies in the 
redshift range 0.4 < z < 0.8 to approximately match 
the redshift distribution of SLACS sources, resulting in 
3369 galaxies. These are shown as small black points or 
hatched blue histograms in Figs. [3HH1 We note that this 
GEMS subsample represents the brighter part of the 
GEMS full HST catalog with surface photometry, owing 
to the shallower COMBO- 17 multiband photometry 
used to determine photometric redshifts (completeness 
limit of F850LP - 23.5 vs. F850LP ~ 24.5 for the 
HST catalog; Barden et al. 2005). 

We computed stellar masses for al l objects i n the 
GEMS sample by applying the recipe bv lBell et all ([2003L 
for a Kroupa IMF) to the rest frame photometry pro- 
vided in the GEM S catalog. As discus sed before, for the 
SLACS sample the I Auger et al . (2009) stellar masses are 
in excelle nt agreement wit h the ones obtained using the 
recipe bv lBell et all ([20031 ). ensuring; that we are able to 
make this comparison. 

4.2. HUDF -selected comparison sample 

For the second comparison sample, we use the galax y 
catalog from the HUDF analyzed bv lCoe et all (|2006D . 
We selected galaxies in a manner similar to the process 
used for GEMS. Galaxies included in our sample are 
those: (1) with g ood G ALFIT fits, (2) that are matched to 
iBeckwith etHI (|2QQ3h ACS objects within 0'/5 (3) with 
successful Bayesian photometric redshifts and (4) which 
are detected in B, V, i', and z' . Our criterion for success- 
ful galfit and redshift fits is that the x 2 be within 4.5a 
of the mean. We also restrict our sample to galaxies in 
the redshift range 0.4 < z < 0.8 to match the redshifts of 
the GEMS and SLACS source samples, for a final tally of 
841 galaxies. These are shown as black points or hatched 
red histograms in Figures BHH1 

We use kcorrect from iBlanton fc Roweisl (|2007[ ) to 
compute rest-frame luminosi ties of the HUDF compari- 
son sample. We then use the lBell et all (|2003l ) algorithm 
to calculate stellar masses, as described in Section I47TT 

4.3. Distribution of observed properties 

We plot the distributions of source Sersic index and 
apparent F814W source magnitude in Figure [3j showing 
both the SLACS sources and the GEMS sample. The 
vertical dotted lines show the magnitudes tested by our 
simulations, demonstrating that we are in a regime where 
source properties can be reliably measured by klens. In 
this figure and throughout this work, the quoted source 
galaxy magnitudes are for the unlensed galaxy, i.e. cor- 
rected for magnification. 

When considering apparent magnitudes, the SLACS 
sources lie squarely between the two regimes sampled 
by the GEMS and HUDF surveys. Most SLACS source 
galaxies fall in the range 22 < F814W < 26, peaking at 
magnitudes fainter than the completeness limit of the 
GEMS survey but brighter than most galaxies in the 

10 The HUDF catalogs used in this work can be 
obtained from the original paper or from D. Coe at 
http : //adcam . pha . j hu . edu/ coe/UDF/. 
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HUDF sample. As mentioned earlier, the drop-off of the 
GEMS catalog at - 23.5 is mostly due to the COMBO- 
17 completeness limit. However, a substantial fraction of 
the SLACS sources are fainter than even the GEMS HST 
completeness limit, suggesting that we are indeed explor- 
ing a different population of intrinsically fainter objects. 
The HUDF survey is much deeper than the SLACS sur- 
vey (144 HST orbits versus one in I band), so it is not 
surprising that these objects are fainter than the SLACS 
sources. The average magnitudes for the SLACS, GEMS 
and HUDF surveys are 24.3, 22.0 and 27.5 respectively. 

In contrast, the distribution of Sersic indices is essen- 
tially the same for the three samples - peaked at n ~ 1 
generally interpreted as dominated by faint disks or com- 
pact galaxies - even though the magnitude ranges probed 
by the three surveys are very different. The median for 
both SLACS and HUDF is 1.1; for GEMS it is 1.3. 

The difference between the SLACS sources and the 
GEMS galaxies is even more pronounced in terms of the 
size distribution, shown in Figure |4] with the same nota- 
tion as Figure [3l The size distribution of GEMS sources 
is much broader and peaks at 0'/78, petering off below 
0'/2. In contrast, Figure |4] highlights the similarity of the 
SLACS and HUDF galaxies. The overall size distribu- 
tions for the two samples is remarkably similar, though 
the HUDF sample extends to much larger sizes, beyond 
the plotted region. This long tail pushes the average 
size up to 0"35 (2.3 kpc) for HUDF, compared to Of! 19 
(1.24 kpc) for the SLACS sources. The median sizes, 
however, are more similar: / /12 for HUDF and / /14 for 
SLACS, corresponding to 0.8 kpc for both samples. 

The complementarity of the SLACS, GEMS and 
HUDF samples is further illustrated by Figure [5] where 
sizes are plotted against apparent magnitude. With 
one exception, SLACS sources have comparable surface 
brightness to that of the GEMS sources, but extend much 
further down in magnitude and size. A few of the bright- 
est SLACS sources have magnitudes comparable to those 
of the GEMS sources, but have sizes at the compact end 
of the distribution. In contrast, the SLACS sources have 
significantly higher surface brightness on average than 
the HUDF galaxies. 

In Figure [5j crosses represent galaxies in the low red- 
shift bin (0.4 < z < 0.6) while squares indicate the higher 
redshift galaxies (0.6 < z < 0.8). The galaxies in all three 
samples are not strongly segregated by redshift, consis- 
tent with a broad distribution in intrinsic luminosity. 

The comparison makes it apparent that the SLACS 
source population is significantly different from the one 
studied by typical HST broad band imaging surveys such 
as GEMS. These sources are most likely excluded by HST 
field surveys because they are too faint and/or too com- 
pact to be detected and identified as galaxies (as op- 
posed perhaps to stars). Furthermore, even if detected, 
they are generally too faint and compact to determine 
their surface brightness profile, size and stellar mass. The 
HUDF survey relies on its extremely large exposure times 
to detect and model such sources. Conversely, SLACS 
identifies them as galaxies and measures their redshifts 
via their emission lines (even though they are detected 
in all three broad bands), and therefore their intrinsic 
faintness and compactness does not represent an obsta- 
cle to detection, while magnification helps in determining 
their structural parameters. In fact, the distribution of 



SLACS sources and sizes is in agreement with a generic 
feature of lensing survey. Owing to magnification bias, 
faint and compact sources tend to dominate the source 
population when their numbe r density de creases sharply 
with luminosity and size (e.g. Treu 2010, and references 
therein). We will return briefly to the lensing selection 
function in the next section. 

4.4. Distribution of rest-frame quantities 

We now explore the distribution in size, luminosity, 
color and stellar mass for the SLACS source galaxies. We 
use two-dimensional plots to examine bivariate distribu- 
tions in the space of well known correlations representing 
physical mechanisms: the size-luminosity relation, the 
color-magnitude diagram, and the size-stellar mass rela- 
tion. Again we use the GEMS and HUDF samples as a 
comparison. 

We begin with the size-luminosity relation in Figure [6] 
We fit a size-luminosity relation of: log 10 r eff (kpc) = 
—0.12 x V — 2.25. However, this result should be in- 
terpreted with a care, keeping in mind the small num- 
ber of objects in our sample and the uncertain selec- 
tion function. A linear extrapolation of the GEMS size- 
luminosity relation shows little offset between the GEMS 
and HUDF galaxies, but the SLACS sources are offset to- 
wards smaller sizes by 0.6 kpc, or ~ 30%. This offset may 
indicate that there is an extended tail of ultracompact 
galaxies of which the SLACS sources are a part, simi- 
l ar to the population of compact em i ssion line galaxies 
(IKoo et all 119951: iPhillips et al.l ll997l iDrozdovskv et al.l 
BOOS IStraughn et aLll2Q09F ~ 

We note that, as expected for a flux limited sample, 
the intrinsically fainter objects are only detected at the 
lowest redshifts (in Figs. [6JHJ crosses represent galaxies 
between 0.4 < z < 0.6 and squares those between 0.6 < 
z < 0.8). Other than this luminosity segregation, no 
evolution is apparent within our dataset. 

The color-magnitude diagram shown in Figure [3 paints 
a similar picture of the properties of the SLACS sources. 
They span the range of blue colors typical of star- forming 
galaxies, and of the GEMS and HUDF samples, even 
with their different intrinsic luminosities. Again, other 
than luminosity segregation, no evolution is apparent. 
Our data (and GEMS) are also consistent with no re- 
lation between color and magnitude, in contrast to the 
sample of HUDF galaxies we have selected for compar- 
ison, which do show bluer galaxies at higher redshift. 
This perhaps surpri sing result is consistent with those of 
iPirzkal et al.l (|2006f ). who looked at emission line galax- 
ies (ELGs) in the HUDF, with -U < M B < -22 and 
rest-frame B-V colors from ~ to ~ 1. As we discuss 
below, the SLACS sources may be best represented as 
ELGs, which could be the cause of this similarity. 

Finally, the offset in size with respect to GEMS and 
HUDF is pronounced even in the size-stellar mass plane 
(Figure [8j), where the SLACS sources seem to define a 
sequence offset to smaller radii by 0.5 kpc, or ~ 30%. 
This indicates that the offset in the size-luminosity re- 
lation is not due to abnormally low mass-to-light ra- 
tios but to an intrinsic compactness of the population. 
For the SLACS sources, we find a size-mass relation of 
logio r eff(kpc) = 0.24 x log 10 M*/M© - 2.20, although we 
again caution the reader as to the small sample size and 
uncertain selection function. 
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Fig. 3. — Characteristics of the SLACS lensed sources. The left panel shows the distribution of Sersic indices for our galaxies (solid unfilled 
histogram), for GEMS sample with 0.4 < z < 0.8 (blue, hatched histogram) and for the HUDF sample with 0.4 < z < 0.8 (red, hatched 
histogram). The right hand panel shows the distribution of source apparent magnitudes, corrected for magnification. To compare with the 
GEMS and HUDF samples we have assumed that F814W=F850LP=F775W (AB), for simplicity. The vertical dashed lines indicate the 
magnitudes in which KLENS was tested by simulations. 
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Fig. 4. — Distribution of source effective radii for our galaxies. The left-hand panel shows sizes in arcseconds; the vertical dotted lines 
indicate the sizes tested by simulations. The right hand panel shows sizes in kpc. The SLACS sources sample is shown as solid histograms; 
the GEMS and HUDF comparison samples are shown as hatched blue and red histograms respectively. 



We note that the stellar masses of the SLACS sources 
are remarkably low for a population at cosmological dis- 
tances, extending below 10 8 Mq. Three systems in par- 
ticular have remarkably low stellar mass: 
J1538+5817 O s =0.531, M*=1O 7O3 M , L y =lO 79 L ), 
J1402+6321 (z s =0.481, M*=1O 6 "M , L y =lO 78 L ), 
J1719+2939 (> s =0.578, M*=1O 686 M , Ly=lO 7 7 L ). 

Their V-band luminosities are just 0.5-0.8 dex higher 
than those of the brighter Milky Way dwarf satellites, 
Sagittarius and Fornax, and conside rably fainter than the 
Large Magellanic Cloud (3 • 1O 9 L Ivan der Mar el et al.l 
As can be seen in Figure [TJ these three sources 



Large 



are being lensed into high magnification partial Einstein 
rings (their total magnifications are the three highest: 
44.3, 31.2 and 40.7 respectively). The J1538+5817 source 
may be a low luminosity satellite of a brighter compan- 
ion, visible just outside the ring. 

4.5. SLACS sources as emission line galaxies 

Although the comparison between our galaxies and 
GEMS or HUDF is useful - it allows us to compare 
the physical properties of our galaxies, including stellar 
mass, to large samples with the same redshift distribu- 
tion - the selection algorithm for SLACS is quite different 
from that of GEMS or HUDF. The SLACS sample was 



Lensed galaxies 



11 



10 



o.i 



0.01 




® a GEMS, HUDF 0.4<z<0.l 
□ SLACS sources z<0.6 
x SLACS sources z>0.6 



J_ 



_L 



1.5 



30 



25 

F814W (AB) 



20 



Fig. 5. — Size-magnitude relation for SLACS source galaxies 
(large points with error bars, with values corrected for mag- 
nification) and all GEMS and HUDF galaxies (small points). 
F814W=F850LP=F775W (AB) has been assumed in the compari- 
son for simplicity. Crosses represent galaxies between 0.4 < z < 0.6 
and squares those between 0.6 < z < 0.8. The diagonal line rep- 
resents constant surface brightness of 24 mag arcsec - 2 in F814W 
(AB), close to the completeness limit of GEMS for a circular source. 
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Fig. 7. — Color magnitude relation for SLACS source galaxies 
(large points with error bars bars) and all GEMS and HUDF galax- 
ies (small points). Crosses represent galaxies between 0.4 < z < 0.6 
and squares those between 0.6 < z < 0.8. We use rest- fram e B and 
V band magnitudes calculated as described in Section RPI For our 
galaxies, we show the rest-frame magnitudes calculated from our 
stellar mass fitting code. 
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Fig. 6. — Size-luminosity relation for SLACS source galaxies 
(large points with error bars bars) and all GEMS and HUDF galax- 
ies (small points). Crosses represent galaxies between 0.4 < z < 0.6 
and squares those between 0.6 < z < 0.8. 



selected ba sed on the detectio n of source galaxy emis- 
sion lines ([Bolton et al.ll200l : GEMS and HUDF are 
both imaging surveys. It is therefore useful to compare 
the SLACS sample with emission line selected galaxies 
(ELGs). SLACS is sensitive to observed fluxes as low as 
-6-10~ 17 erg s _1 cm- 2 (jBolton et all 120041 ; IBobler et all 



120081 ). Taking into account the effects of magnification 
the equivalent depth is typically ~7-10 _18 erg s _1 cm -2 , 
but can be higher by a factor of a few for the most ex- 
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Fig. 8. — Size-mass relation for SLACS source galaxies (large 
points with error bars bars) and all GEMS and HUDF galaxies 
(small points). Crosses represent galaxies between 0.4 < z < 0.6 
and squares those between 0.6 < z < 0.8. Stellar masses for GEMS 
and HUDF galaxies have been computed using the Bell et al. 2003 
recipes. All stellar masses have been converted to a Kroupa IMF 
normalization. 

treme objects. 

Emissio n line galaxies identified by HST grism sur- 
veys (e.g. iDrozdovskv et al.l 120051 : IStraughn et all 120091 ) 
reach comparable depths, and provide a good bench- 
mark. Those samples span similar broad band magni- 
tude ranges to the SLACS sources and have compara- 
ble sizes. For ex ample, the size dis t ributi on of galaxies 
in the sample of IDrozdovskv et al.l ([200 5) is limited to 
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small sizes (r e ff < Of! 5) and peaks at r e ff ~ O'/l — 0'/2. 
Magnitudes peak at F814W~24, reaching magnitudes as 
faint as 27. These are quite similar distributions to those 
for the SLACS galaxies. 

5. DISCUSSION AND CONCLUSIONS 



Following an initial foray by Marshal l et al.l ([20071 ), 
we have extended the study of the size-mass relation 
of galaxies to lower masses and smaller sizes by mak- 
ing use of the factor of ~ 10 magnification provided by 
the gravitational lenses of the SLACS survey. In order 
to model this unprecedented sample of 46 gravitational 
lens systems in three bands (Vhst, ^814? #160 )> we have 
developed and tested klens, a new code for fast lens 
modeling optimized for this work. Extensive testing and 
comparison with standard codes like galfit show that 
we are able to accurately measure the properties of small, 
faint lensed galaxies: for a mock lensed galaxy of 25th 
/8i4-band magnitude, we are able to recover a Of! 05 size. 

The main result of this study is that by exploiting the 
gravitational lensing effect we have been able to deter- 
mine sizes and stellar masses of a sample of very faint 
and compact galaxies, which are typically not studied 
by imaging surveys and are beyond the reach of normal 
spectroscopic follow-up. We can derive accurate sizes 
and stellar masses for SLACS sources that are typically 
1-2 magnitude fainter and 5 times smaller than those of 
a galaxies selected for structural analysis by typical HST 
imaging survey (GEMS). Although the lensed sources are 
not as faint as galaxies from extremely deep imaging sur- 
veys (HUDF), we are able to reach /gi4 magnitudes of 27 
with limited exposure times. For comparison, the HUDF 
survey is comprised of 400 HST orbits, with 144 in I- 
band, while the SLACS survey is just one orbit in each 
band. Furthermore, we identify galaxies that are typi- 
cally ~ 30% smaller in effective radius than those iden- 
tified in the HUDF at comparable luminosity or stellar 
mass. 

The SLACS sources are similar to emission line 
selected-galaxies (e.g 



ll997HDrozdovskv et all 



Koo et all Il995t Phillips et al.l 



20051 : IStraughn et alJ l2009). and 



as such may represent the building blocks of present day 
dwarf-spheroidals or perhaps even bulges of galaxies (see 
lHammer et al.ll2001f ). However, thanks to the lensing ef- 
fect we can image them at high angular resolution and 
derive accurate sizes and structural parameters. 

The existence of this population highlights the impor- 
tance of diverse selection criteria for a complete census 
of the galaxy population. Lensing imposes a completely 
different selection function than for imaging surveys: in 
particular, the SLACS sources were selected by their 
strong, multiple emission lines, which were observed to 
stand out against the background of a n early-type (lens) 
galaxy spectrum ([Bolton et al.ll2004l ). In contrast, the 
GEMS and HUDF sources were selected on the basis of 
their broad band flux and morphology, using no informa- 
tion as to their brightness in emission lines. Similarly, it 
is likely that spectroscopic lens surveys could preferen- 
tially select compact galaxies, owing to the larger mag- 
nification for compact sources. Modeling the selection 
fu nction in detail (for example along the lines suggested 
by iDobler et al.l [2008) is essential to quantify the lumi- 
nosity function of this population and its abundance rel- 
ative to the population of broad-band selected galaxies, 



and that of emission line selected galaxies. This model- 
ing work goes beyond the scope of this paper but will be 
pursued in the future with the goal of determining the 
shape an intrinsic scatter of the size-mass relation taking 
into account the lensing-selected population. 

A brief summary of the main results of this paper fol- 
lows: 

• Our best-fit lens mass parameters are in good 
agreement with those found in Paper V and Pa- 
per IX in the Jgu band. To mitigate against the 
scatter introduced by lens light subtraction in the 
other, lower S/N filters, we fixed the mass models 
to those fitted in the highest S/N band. 

• The inferred sources appear to be similar in struc- 
ture to the g alaxie s in the GEMS catalog of 
Bar den et all (|2005D and the HUDF catalog of 
Coe et all ( 2006f h in the sense that the distribu- 



tions of Sersic indices are the same between the 
three samples. However, the SLACS sources have 
a significantly different distribution of magnitudes. 
They are fainter than GEMS, with the peaks of 
their differential number counts separated by 2 
magnitudes, and brighter than HUDF by about 4 
magnitudes. (However, SLACS reaches these mag- 
nitudes with only one HST orbit in I-band, while 
HUDF uses 144). Aside from this shift, the color- 
magnitude diagrams of the three samples are con- 
sistent with each other. 

• The SLACS sources are also significantly smaller 
than the GEMS galaxies and have sizes similar to 
the HUDF galaxies. The SLACS sources are made 
measurable via the lens magnification, while the 
extremely long exposure time is what allows study 
of the HUDF galaxies. The size distributions of 
the SLACS sources and the HUDF galaxies peak 
at approximately 0.1 arcsec (~ 0.8 kpc); such small, 
faint objects were likely below the COMBO- 17 flux 
limit used in constructing the GEMS sample, or 
were unresolved in the shallow HST images. 

• The closest analog in size and magnitude to the 
SLACS sources are the faint emission line galaxies 
identified by blind grism surveys with HST. 

• Combining the reconstructed (model) magnitudes 
from the three HST bands, we infer stellar masses 
for each SLACS source and plot the size-mass re- 
lation at 0.4 < z < 0.8. We then compare to the 
GEMS and HUDF results, finding that the SLACS 
sources are offset from the GEMS and HUDF size- 
mass and size-luminosity relations towards smaller 
effective radii. 

• Some of our measured objects are very low stellar 
mass indeed, ~ 1O 7 M , comparable to the largest 
dwarf satellites of the Milky Way and fainter than 
the Large Magellanic Cloud. 

This work is a further step towards extending the study 
of intermediate and high redshift galaxies to smaller and 
fainter galaxies, demonstrating that gravitational lens- 
ing can help us probe this regime. A larger sample of 
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lenses will improve the statistical significance; the anal- 
ysis code we have developed can perform lens modeling 
simply, quickly, and robustly and will be useful when 
even large samples of lenses are available. Before this, a 
fuller understanding of the lensed source selection func- 
tion will allow us to make a quantitative analysis of the 
size-mass relation of dwarf galaxies at intermediate red- 
shift and how this relation evolves with cosmic time. 
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APPENDIX 

MEASURING SOURCE GALAXIES WITH klens 

We aim to measure the magnitude and effective radius of the source galaxies of the gravitational lenses in our sample. 
In this section we describe a new code, klens, designed for fast and robust lens modeling of large number of images 
(for this paper alone we need to model 138 indepen dent images, i.e. 3 bands for each of the 46 lenses). It is similar in 
spirit to the software used by Marsha ll et al.l (j2007[ ). in that the models for the lens mass and source surface brightness 
distributions are the same. However, we do not explore the posterior probability distribution (PDF) for the model 
parameters fully, but instead simply search for the peak of the distribution. This approach saves a considerable amount 
of CPU time (up to two orders of magnitude), since the computation of the image pixel data likelihood function is 
relatively expensive - but it gives us only the covariance matrix in the Laplace approximation as opposed to full 
statistical uncertainties on the parameters. We note however, that the error budget in this kind of analysis is totally 
dominated by systematic errors due to factors such a s lens light subtraction, model simplicity and PSF approximation 
(see the following sections, and Marshall et al. 2007). Therefore we conclude that, with current CPU limitations, it is 
more effective to combine minimizatio n with Laplace appro ximation of the errors with a detailed study of one or more 
objects to estimate systematic errors (M arshall et al.ll2007[ ). 

After describing the lens and source models in a bit more detail, we outline our implementation of the prior PDFs 
and the optimization of the posterior, before demonstrating the code's performance on three test datasets. 

Models 

The gravitational lense s in our sample are kn own to be well- approximated by simple, singular isothermal ellipsoid 
(SIE) mass distributions (K oopmans et al.ll2QQ6h . We therefore use an SIE model to describe the mass distribution of 
the lens galaxy. This model has five para meters: the velocity d ispersion <tsie, elliptical axis ratio q = b/a, orientation 
angle and position (x&,y&). (See e.g. iKormann et aHfTODl. for details). For the source galaxy surface brightness, 
we use an elliptically symmetric distribution with Sersic profile. The source model has seven parameters: the position 
(x s ,y s ), orientation angle # s , inclination angle i (more often referred to be cosz, which equals the apparent ellipse axis 
ratio), effective radius r e fr, AB apparent magnitude m, and Sersic index n. A Sersic ind ex of n = 0.5 cor responds to a 
Gaussian, n = 1 to an exponential disk, and n = 4 to a de Vaucouleurs profile (see e.g. iPeng et "aTll2002L for details). 

In Paper V we used multiple Sersic profile components to describe the source galaxy. Because our goal here is to 
measure the physical properties of the source such as size, magnitude and mass, and compare them with corresponding 
measurements of non-lensed galaxies in the literature, we choose instead to describe the source as a single object. As 
a result, we do not expect to be able to fit all our sources to the same level of detail as we did in Paper V. One option 
would be to remove all lenses with complex sources from our sample; however this would bias us towards smaller and 
simpler sources, which is undesirable. This same single-component simplification was used in, for example, the GEMS 
and SDSS size-mass relation studies; restricting ourselves to single-component sources allows us to to make direct 
comparisons with these nonlensing surveys. 

Asserting that the a ssumption of an SIE m odel for the lens galaxy mass profile is the most significant source of 
systematic uncertainty, [Marshall et al.l (|2007| ) estimated the systematic errors in the source size and brightness to be 
approximately 12% (0.07 kpc for a 0.6 kpc source at z s =0.6) and 0.26 magnitudes respectively. However, we note that 
these errors will be common to all observations of a given lens. 

Priors 

To ease the exploration of the model parameter space, we work with hidden parameters constrained to vary uniformly 
on a hyper cube, such that their typical values are around one. These hidden parameters are then transformed to the 
physical parameters described in the previous section before the predicted image is generated. This transformation 
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TABLE 3 

Model transformations and hard limits 





"m pfpr 

_L CX± CX111C- IjC-± 






Limits 




Axis ratio q = b/a 


Uniform 


Linear 


[0,1] 


Lens mass profile 


Velocity dispersion csie (km s — 1 ) 


Uniform 


Linear 


[150, 400] 


Orientation angle 6& (rad) 


Uniform 


Linear 


[0,7T] 




Position Xd,2/d (arcsec) 


Gaussian 


Error function 


[-1-1] 




Position angle (rad) 


Uniform 


Linear 


[0,7T] 




cos(i) = q 


Uniform 


Linear 


[0,1] 


Source surface brightness profile 


Effective radius r e Q (arcsec) 


Uniform 


Linear 


[io- 3 ,ioo] 




Magnitude m 


Power law, index k 


Power law , index k 


[10, 30] 




Sersic index n 


Uniform 


Linear 


[0.2, 10] 



need not be linear: in fact, it is used to encode our prior knowledge about the model parameters. For example, a 
uniform prior PDF is implemented as a linear transformation, while exponentiating a hidden parameter corresponds 
to assigning a scale- free ("Jeffreys", 1/x) prior. Truncation of the prior PDFs is achieved where necessary by setting 
the log likelihood to be very large and negative outside the specified limits. The parameters of our model are listed 
in Table [3j along with the prior PDFs assigned to them, unit transformation used, and the range over which they are 
allowed to vary. 

Posterior evaluation and optimization 

At each point in parameter space visited we compute the pixel value likelihood function. We assume Gaussian 
errors on the pixel values, given by the weight (inverse variance) map produced during data reduction (Paper IX). The 
log likelihood is then just minus half the image x 2 . The predicted image, required for this statistic, is computed by 
mapping each pixel's position back to the source plane via the lens equation, and looking up the value of the model 
source surface brightness. We over-sample the image p lane to mak e sure we resolve the source. The lensed images are 
then convolved with a model PSF; we use Tiny Tim (Krist 2003) to generate this PSF. Note that we do not model 
lens galaxy light, only the source surface brightness. We therefore assume perfection subtraction of lens galaxy light, 
an assumption which incurs the systematic error discussed at the end of Section [2] 

Since the prior PDF is implemented by parameter transformation, optimization of the posterior is equivalent to 
maximi zing the log likeli hood with respect to the hidden parameters. For this optimization we use the IDL routine 
mpfit (jMarkwar dt 2009), which is a generalized implementation of the Levenberg-Marquardt nonlinear least-squares 
algorithm. 

The likelihood function is typically quite sharply peaked, with multiple local maxima to be avoided. We find that a 
three- step process is required for the global peak to be successfully located. First, we make an initial estimate of the 
source galaxy parameters using a best-guess mass model: the image pixels are traced through the mass model using 
the lens equation and flux is accumulated on a grid in the source plane. Where multiple image pixels map to the same 
source plane pixel, we take their simple average. We then take the initial source galaxy position to be the first moment 
of this resulting light distribution, while summing all the source plane pixel values provides an initial approximation to 
the source magnitude (this is usually an over-estimate). We gain no knowledge of the source inclination, size, position 
angle or Sersic index in this first stage; the source position and magnitude are the most important for finding a good 
fit to the data, in the sense that small deviations in these parameters can lead to very low likelihood models. 

Next, we refine our initial estimates by running mpfit, using our initial estimates of the source position and mag- 
nitude as its starting point (we initialize the other parameters at typical values). This preliminary run is performed 
without image plane subsampling, for speed. We also use a deflated weight image (the square root of the weight) in this 
second stage; this smooths the likelihood function, making it easier for the optimizer to locate the global maximum. 
Finally, in the third step we re-start mpfit at the position found in step two to get final maximum posterior model 
parameters for the mass and source. For this final run, we use an image plane over-sampling of 4 (along each axis), 
and the normal weight image. Convolution with the PSF is done on the over-sampled grid; we then use neighborhood 
averaging to return to the original size. 

Testing klens 

In order to verify that the klens code works properly, we performed tests on (1) simulated non- lensed galaxies, (2) 
real non-lensed galaxies and (3) simulated lenses. All simulated data was created using HST/ ACS igu-band (F814W) 
image parameters and PSFs, while the real galaxies were selected from HST/ ACS /gu-band images. 

Testing klens on simulated non-lensed galaxies 

We first test klens on simulated non-lensed galaxies. Each mock galaxy is taken to be a single elliptically symmetric 
Sersic profile component, to allow us to investigate straightforward parameter recovery by the code. The sizes, 
magnitudes and Sersic indices of these galaxies were varied over a 3 by 3 by 3 grid: we chose sizes from (O'/IS, 0'/3, 0'/6), 
magnitudes from (20.5,21.5,22.5) and Sersic indices from (0.5,1,2). All position angles were set to S = 1.75 and 
inclinations at cos(z) = q = 0.6. The simulated galaxies were placed at the center of 42 x 42 pixel images. Magnified 
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TABLE 4 

Mean parameter difference (inferred - true) and associated rms scatter for simulated non-lensed galaxies 



Parameter KLENSmodels GALFlTmodels 



(Am) 


(118 ±1) x 10- 4 


(-76.0 ±0.6) x 10- 4 


Cm 


0.027 


0.016 


(Ar eff > 


(-730" ± 7) x 10- 5 


(546" ± 4) x 10- 5 


<j r 


0.019 


0.012 


(An) 


(-350 ±3) x 10- 4 


(159 ± 2) x 10- 4 




0.079 


0.041 




21 22 
magnitude 



23 



£0.4 



0.2 




21 22 
magnitude 



23 



Fig. 9. — Testing klens and galfit on simulated non-lensed galaxies. Best-fit effective radii and magnitudes for klens (left panel) and 
GALFIT (right panel), binned by true size and magnitude: the dashed vertical and horizontal lines indicate these bins. Errors bars indicate 
the r.m.s. scatter in each bin (the errors on the mean are ~ 5.5 times smaller and are not shown). In the largest size bins, klens tends to 
underestimate size and overestimate magnitude, while GALFIT tends to do the opposite. 

by a typical lens, this size is approximately equivalent to the 122 x 122 pixel cutouts used for the lenses. We generated 
ten independent noise realizations for each galaxy, giving a total of 270 simulated observations. The whole sample 
were then modeled with both klens (using the priors given in Table [3]) and galfit. 

The results of the parameter recovery are given in Table |4j which lists the mean (inferred - truth) difference, and 
associated r.m.s. scatter values. We look at the results more closely in Figure [9j Overall klens and galfit give 
very similar results, although the rms scatter is larger with klens. For both modeling programs, the largest galaxies 
contribute most of the bias, which is in the sense that klens tends to slightly underestimate sizes, while galfit tends 
to slightly overestimate them, even though the residual images are consistent with simulated noise. This is because the 
two programs are sampling different parts of the covariance between size, magnitude and Sersic index. To illustrate 
this point, we plot the best-fit Sersic indices against the best-fit sizes in Figure [lOj again grouped by true parameter 
value. 

Test on real non-lensed galaxies 

Next, we tested klens on real non-lensed galaxies. We selected a sample of 981 galaxies, detected in a subset of the 
SLACS fields, to match the nine magnitude and size bins used in the previous test: the sample consists of all objects in 
those fields that have SExTRACTORmagnitudes 20 < m < 23 and O'/l < r < 0'/65 that are not stars or image defects. 
We then ran klens on each test image, with the priors given in Table [3] we also ran galfit on each dataset. 

For this sample we do not know the true values for each object's size, magnitude and Sersic index; we can only 
compare galfit and klens results. When comparing Sersic indices, we noted significa nt sys tematic differences 
between the two programs, especially at n > 1.5. This scatter is, as discussed in Section [A.4.1[ a result of klens 
and galfit sampling different parts of the covariance between Sersic index, size and magnitude. We compare the 
two codes' best-fit magnitudes in Figure [TTJ In this plot we highlight those galaxies for which best-fit Sersic indices 
differ by less than 20%; 408 objects meet this requirement. The magnitudes and sizes for these galaxies were found 
to to be in reasonable agreement: (Am) = 0.147 ± 0.008 with r.m.s. scatter 0.170. A similar analysis for size gives 
(Ar e ff ) = —0.174 ± 0.009 with r.m.s. scatter 0.175. As we decrease the allowable difference in Sersic index, the bias 
and scatter both decrease. 

galfit occasionally returns very low magnitudes relative to the SExtractor values, and that these are associated 
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Fig. 10. — Testing klens and galfit on simulated non-lensed galaxies. We plot best-fit Sersic index against the ratio of the best-fit to 
true size, in order to illustrate the covariance between Sersic index, size and magnitude. Galaxies are again binned by true Sersic index 
and size, shown by the dashed lines. Red triangles indicate GALFIT fits and blue squares indicate klens fits. The error bars indicate the 
mean value and scatter of galaxies in the faintest bin (the errors on the mean are ~ 3 times smaller and are not shown). The covariance - 
and the fact that klens and galfit sample disparate parts of it - is most evident in the n = 2, r = 0'/6 bin. 

with Sersic indices very different from the klens values. We note that, for example. iBarden et al.l (j2QQ5h reject objects 
with galfit magnitudes differing from their SExtractor magnitudes by more than 0.6 magnitudes. 

Test on a simulated lens sample 

Finally, we tested the ability of klens to recover the structural parameters of lensed source galaxies. We can do this 
to a limited extent by comparing our klens mass models to the original SLACS le ns mod els (Section 13. 2p . However, 
for the source parameters we need to use simulated galaxies, as we did in Section IA.4.11 in the non-lensed case. We 
generated a sample of lensed galaxies to be representative of the possible lensed image configurations and sources. 
The exact properties of the lens itself are less important: we used the same mass profile (<tsie=270 km s _1 , q=0.77 
and #d=1.75 rad, typical parameters for a SLACS lens) for each synthetic system. We also fixed the source galaxy 
orientation at S =1.75 rad and its inclination at cos(i) = q = 0.6, again typical values. 

The source position, effective radius, magnitude and Sersic index were varied, to give 108 possible lenses. We chose 
four source positions such that the four main lens morphologies (double, caustic, cusp and quad) were represented. 
The source effective radius was chosen from (0'/05, O'/l, 0'/2), the magnitude from (23, 24, 25) and the Sersic index from 
(0.5, 1.0,2.0). Taking into account magnification due to lensing, these sizes and magnitudes are roughly equivalent to 
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Fig. 11. — Testing klens against galfit for real galaxies. Shown are the best-fit magnitudes from galfit, on the x-axis, and from klens, 
on the y- axis. Large filled blue squares are those points for which the percent difference in Sersic index is less than 20%. The red dashed 
lines indicate the upper and lower limits on the SExtractor magnitudes for objects in the sample. 



those used in the non- lensed galaxy tests. For each mock system we simulated pixelated images, again with no lens 
galaxy light. We generated ten noise realizations for each lens, to give a total of 1080 mock observations. We then ran 
klens on each mock dataset, with the priors given in Table [3] 

Parameters for which the same true value was used in all simulations (i.e. the lens parameters, and the source 
inclination and orientation angle) are well- modeled, with only small biases well below the r.m.s. scatter (see Table [5]). 
We divide the sample into bins to consider the inferences of the remaining parameters, their means and scatters across 
the sample; these results are presented in Table [6l In Figure [12] we look at the performance of klens in recovering 
the input effective radius and magnitude, binned by their true values. 

We find that the inferred parameters are consistent with the input values, albeit with larger scatter than was seen 
in the non- lensed simulation case. This is partly due to our use of fainter galaxies, which are harder to measure when 
the image configuration is of lower magnification. (We note that the scatter is smaller for high surface brightness 
galaxies, indicating that - as might be expected - it is easier to model such galaxies.) It may also be due to the way 
the information in the data is being used to simultaneously constrain the lens model - to some extent changes in the 
source parameters can be balanced by changes in the lens model. 
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TABLE 5 

KLENS PARAMETER RECOVERY IN A SIMULATED LENS SAMPLE COMMON LENS AND SOURCE PARAMETERS 



Parameter 


True value 


Inferred value 


r.m.s. scatter 


q 


0.77 


0.76044 ± 0.00006 


0.063 


asm (km s -1 ) 


270 


269.961 ± 0.007 


7.294 


d (rad) 


1.75 


1.7368 ±0.0002 


0.221 


S (rad) 


1.75 


1.7288 ±0.0003 


0.293 


cos(z) 


0.6 


0.58685 ± 0.00004 


0.044 



TABLE 6 

KLENS PARAMETER RECOVERY IN A SIMULATED LENS SAMPLE VARIABLE SOURCE PARAMETERS 



Parameter 


True value 


Inferred value 


r.m.s. scatter 




0.05" 


0.0543" ± 0.0006 


0.011" 


effective radius 


0.1" 


0.1023" ±0.0009 


0.016" 




0.2" 


0.195" ±0.001 


0.022" 




23 


23.01 ±0.02 


0.107 


magnitude 


24 


24.00 ±0.01 


0.062 




25 


24.95 ± 0.02 


0.265 




0.5 


0.62 ±0.02 


0.354 


Sersic index 


1.0 


1.05 ±0.01 


0.247 




2.0 


1.98 ±0.02 


0.291 



Summary 

We have tested KLENS on three different types of data. Our results on simulated, non-lensed simple Sersic profile 
galaxies show that KLENS gives almost identical results to galfit; there is a small (< 5% for 0'/6 size objects) systematic 
bias towards larger inferred sizes and Sersic indices with galfit, and a bias of comparable magnitude towards smaller 
inferred sizes and Sersic indices with KLENS, as they sample different parts of the covariance. In real, non-lensed 
galaxies, this covariance results in significant, systematic differences between KLENS and galfit inferred parameters, 
with galfit suffering more catastrophic magnitude errors - these have simply been removed in the literature. The 
difference between the codes is less pronounced if only models with matching Sersic indices are considered. Finally, 
testing KLENS on simulated lens systems where the source is a simple Sersic profile galaxy shows that, for typical 
expected sizes and magnitudes, KLENS is able to provide accurate measurements of these quantities, with its precision 
improving with source surface brightness. From this last test (Figure H2j). we estimate that the statistical uncertainties 
on the source r e ff and brightness are 5% and 0.1 magnitudes respectively. Adding these to the other systematic errors 
noted above in quadrature, we find overall systematic errors of 13% on r e ff, and 0.3 magnitudes photometric error 
(although 0.26 mag in quadrature of the latter - due to assuming an SIE for lens mass model - is common to all filters 
studied). 
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Fig. 12. — Testing klens source parameter recovery on simulated lensed galaxies. We plot KLENS-inferred model magnitude (x-axis) and 
effective radii (y-axis) for a sample of 1080 simulated lenses. The faint dashed lines again indicate the true parameter values (the 9 possible 
combinations of true size and magnitude occur at their intersections). The data is binned by true size and magnitude; the nine points show 
the average values for each bin. The error bars again indicate the r.m.s. scatter of the inferred parameters (the errors on the mean are 11 
times smaller and are not shown). 
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